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Abstract This article is concerned with a stochastic multi-patch model in which each local 
population is subject to a strong Allee effect. The model is obtained by using the framework 
of interacting particle systems to extend a stochastic two-patch model that has been recently 
introduced by Kang and the author. The main objective is to understand the effect of the ge- 
ometry of the network of interactions, which represents potential migrations between patches, 
on the long-term behavior of the metapopulation. In the limit as the number of patches tends 
to infinity, there is a critical value for the Allee threshold below which the metapopulation 
expands and above which the metapopulation goes extinct. Spatial simulations on large reg- 
ular graphs suggest that this critical value strongly depends on the initial distribution when 
the degree of the network is large whereas the critical value does not depend on the initial 
distribution when the degree is small. Looking at the system starting with a single occupied 
patch on the complete graph and on the ring, we prove analytical results that support this 
conjecture. From an ecological perspective, these results indicate that, upon arrival of an alien 
species subject to a strong Allee effect to a new area, though dispersal is necessary for its 
expansion, strong long range dispersal drives the population toward extinction. 



1. Introduction 

To understand the role of dispersal in populations subject to a strong Allee effect, Kang and 
the author recently introduced deterministic and stochastic two-patch models [1]. In population 
dynamics, the term Allee effect refers to a certain process that leads to decreasing net population 
growth with decreasing density. This monotone relationship may induce the existence of a so- 
called Allee threshold below which populations are driven toward extinction, a phenomenon which 
is referred to as strong Allee effect. In this paper, we continue the analysis initiated in [4J and 
use the framework of interacting particle systems to extend the stochastic two-patch model to a 
more general multi-patch model. Thinking of the set of patches as the vertex set of a graph in 
which each edge indicates potential dyadic interactions between two patches, the main objective 
is to understand how the geometry of the network affects the survival probability of the global 
metapopulation. Let G = {V, E) be a graph representing the network of interactions. The system 
is a continuous-time Markov chain whose state at time i is a function 



rjt :V ^ [0, 1] with 'rjt{x) = population density at vertex x. 

Having an Allee threshold 9 € (0, 1) and a migration factor ^ € (0, 1/2], the evolution consists of 
the following two elementary events occurring in continuous time: 

• Mixing events - Each edge becomes active at rate one, which results in a fraction ^ of the 
population at each of the two interacting vertices to move to the other vertex. 



'Research supported in part by NSF Grant DMS-10-05282. 
AMS 2000 subject classifications: Primary 60K35 

Keywords and phrases: Interacting particle system, Allee effect, dispersal, metapopulation. 



1 



2 



N. Lanchier 



• Local events - Each vertex becomes active at rate one, which results in the population density 
at that vertex to jump from below 9 to state or from above 9 to state 1. If the density is equal 
to the Allee threshold then it jumps to either state or state 1 with probability 1/2. 

The inclusion of mixing events indicates that individuals can move from patch to patch through the 
edges of the graph. The strength of dispersal is thus modeled by the mean degree distribution of 
the network of interactions. Local events model the presence of a strong Allee effect in each patch: 
local populations below the Allee threshold are driven toward extinction whereas local populations 
above the Allee threshold expand, and we think of state 1 as the normalized density of a local 
population at carrying capacity. Formally, the Markov generator is given by 

ix,y)£E 

+ (l{r?(x) >9} + (1/2) = 9}) [fiat v) " fiv)] (1) 

+ (Hvix) <0} + (1/2) iMx) = 9}) [f{a- 77) - /(r?)] 
xev 

where (cr+ r/)(x) = 1 and {a~ rj){x) = 0, and where 

(cT^,y ry)(zi) = 7]{zi) + ^ {r]{z2) - -qizi)) whenever {zi, Z2} = {x,y} 

while the state at all other vertices is unchanged. The stochastic two-patch model we introduced 
and studied in [1] is simply the process ([1]) when the network of interactions consists of two vertices 
connected by a single edge. The main objective was to answer the following question: starting with 
one empty patch and one patch at carrying capacity, does the inclusion of mixing events lead to a 
global extinction, i.e., both patches in state eventually, or to a global expansion, i.e., both patches 
in state 1 eventually? Theorem 8 in ^ gives the following answer: when the migration factor is 
small, the probability of global extinction, respectively, global expansion, is close to one when the 
Allee threshold is larger than one half, respectively, smaller than one half. As the migration factor 
increases, the limit becomes less predictable. This result suggests that one half is a critical value 
for the Allee threshold. However, we literally interpreted this one half as one patch initially in 
state 1 divided by two patches, and also conjectured that, for systems in which all the patches 
are connected, i.e., the network of interactions is a complete graph, the critical value for the Allee 
threshold is equal to the initial fraction of patches in state 1 in the limit as the size of the system 
tends to infinity. This conjecture is supported by the last diagram of Figured] which shows numerical 
results for the stochastic process ([T|) on four regular graphs with different degrees. In contrast, the 
first diagram suggests that the critical value for the Allee threshold is again equal to one half 
for the process on the ring in the limit as the size of the system tends to infinity. This indicates 
that the critical value for the Allee threshold strongly depends on the initial configuration in the 
presence of strong long range dispersal whereas the initial configuration is essentially unimportant 
in the presence of weak short range dispersal. The main objective of this paper is to explain at 
least qualitatively this difference between the process on the complete graph, which models strong 
dispersal, and the process on the ring, which models weak dispersal. 
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Figure 1. Simulation results for the process on the torus Z/200Z in which each vertex is connected to its d nearest 
neighbors. Each picture shows the density of vertices below/above the Allee threshold after a large number of updates 
as a function of the Allee threshold (200 values) and the initial fraction of vertices in state 1 (200 values). The color 
for each of the 200 x 200 parameter values is computed from the average of 100 independent realizations, and the 
color code is black for all below the Allee threshold and white for all above. In all the simulations, /i = 0.2. 



2. Main results 

We first assume that the process starts from the product measure in which each vertex is in 
state 1 with probability p and in state otherwise. Following the terminology of [1], we call global 
extinction the event that the process converges to the "all 0" configuration and global expansion 
the event that it converges to the "all 1" configuration. Note that, on finite graphs, the process 
converges to one of these two absorbing states. In particular, the probability of global extinction 
and the probability of global expansion sum up to one, thus showing that the long-term behavior 
is completely characterized by the probability of global expansion 

Pg{0, /i, p) := P{r]t = l for some t), 

which depends on the network of interactions, the Allee threshold, the migration factor, and the 
initial density of occupied patches. The simulation results of Figure [T] suggest some monotonicity of 
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the probability of expansion with respect to the Allee threshold and the initial density of occupied 
patches, as well as a certain symmetry between the probability of expansion and the probability of 
extinction. These results follow directly from standard coupling arguments for interacting particle 
systems that we briefly describe without detailed proof. 

Monotonicity with respect to ^ - Two processes on the same graph and starting from the 
same initial configuration but with different Allee thresholds can be coupled in such a way that 
the process with the smaller Allee threshold dominates the other process, which implies that the 
probability of global expansion pg[9,ijl,p) is nonincreasing with respect to 9. 

Monotonicity with respect to p - Two processes on the same graph, with the same Allee 
threshold, and with the same migration factor can be coupled in such a way that if one process 
dominates the other one at time then the domination remains true at all times. This implies that 
the probability of global expansion pg{0,^,p) is nondecreasing with respect to p. 

Symmetry - The process with Allee threshold 9 can be coupled with the process on the same 
graph and with the same migration factor but with Allee threshold 1 — 9 m. such a way that if at 
any vertex the initial population density for one process equals one minus the initial population 
density for the other process then this remains true at all times. This implies that 



and explains the symmetry in the four simulation pictures of Figure [TJ 

The process starting with a single occupied patch - We now return to the main objective 
of this paper, which is to understand the effect of the geometry of the network on the invadability 
of species subject to a strong Allee effect. This aspect is mathematically more difficult to under- 
stand because two processes on different graphs cannot be coupled in such a way that one process 
dominates the other one. Our analysis focuses on the extreme cases of the ring and the complete 
graph corresponding to the first and last diagrams of Figure [D Let 



where the + and — superscripts allude respectively to the fact that the complete graph is the 
connected regular graph with the largest degree while the ring is the connected regular graph with 
the smallest degree. Recall that the simulation results of Figure [D suggest that 



pg{o,p,p) 



l-pG{l-9,fi,l-p) 



p%{9,n,p) 



the probability of global expansion for the process 
on the complete graph with vertices 

the probability of global expansion for the process 
on the ring with vertices 



(2) 



limAT^oo p%{9,ii,p) 



when 9 > p 

1 when 9 < p, 



(3) 



whereas for the process on the ring starting with p € (0, 1) 



limAT^oo Pj^{9,p,p) 



when 9 > 1/2 

1 when 61 < 1/2. 



(4) 



Following [1] whose main objective was to understand whether an alien species established in 
one patch can expand in space, we assume from now on that the process starts with a single 
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patch in state 1 and all the other patches in state 0. In particular, we drop the parameter p in the 
probabilities ([2])- The following two theorems give qualitative differences between the system on the 
complete graph and the system on the ring starting with a single vertex in state 1, which supports 
the conjectures ([3D and (jH). More precisely, our first theorem indicates that, even when the Allee 
threshold is very small, the system on the complete graph is driven toward global extinction with 
high probability when the number of vertices is large. 

Theorem 1 - Assume that 6*, /i > 0. Then, p'^{9, p,) ^ as N ^ oo. 

In contrast, when the Allee threshold is small enough, the system on the ring expands globally with 
a positive probability that does not depend on the number of vertices. 

Theorem 2 - Assume that 9 < p^ {1 - ;u)^^^°. Then, Mn Pn^^^iA > 0- 

The mysterious assumption in the previous theorem follows from a series of bounds of certain 
probabilities that are estimated based on geometric arguments and are not optimal. Some of these 
estimates appear in our proof and the other ones in the calculation of an upper bound for the 
critical value of one dependent oriented site percolation in, e.g., [1]. Even though the assumption 
of the theorem is far from being optimal, it gives at least an explicit lower bound for the critical 
value of the Allee threshold for the process on the ring. More importantly, the combination of 
both theorems show the following qualitative difference: for some values of the Allee threshold, the 
probability of global expansion is bounded from below for the process on the ring but vanishes to 
zero for the process on the complete graph as the number of vertices increases. This supports at 
least qualitatively the contrast between ([3]) and ([!]). From an ecological perspective, this indicates 
that, upon arrival of an alien species to a new area, though dispersal is necessary for its expansion, 
the best strategy is to first only disperse to nearby patches, and then progressively increase the 
strength of its dispersal as the fraction of patches at carrying capacity increases. 

3. Preliminary results 

This section gives some definitions and simple results that will be used repeatedly in the proof 
of both theorems. Throughout this paper, we think of the process ([T]) as being constructed from 
a Harris' graphical representation [3]. Each edge of the graph is equipped with a Poisson process 
with intensity one while each vertex is equipped with a Poisson process with intensity one and a 
sequence of independent Bernoulli random variables with parameter 1/2. We write 

• Tn{x,y) := nth arrival time of the Poisson process attached to edge {x,y) € E, 

• Un{x) := nth arrival time of the Poisson process attached to vertex x (zV, 

• Bn{x) := nth member of the Bernoulli sequence attached to vertex x . 

All these Poisson processes and Bernoulli random variables are independent and form together a 
percolation structure from which the multi-patch model ([T]) can be constructed. 

• Mixing events - At time t := Tn{x,y), we draw a double arrow along the corresponding 
edge to indicate the occurrence of the following mixing event: 



r]t{x) = {ax,yr]t-)(.x) and r]t{y) = {ax,yVt^){y). 
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• Local events - At time t := Un{x), we put a dot at vertex x to indicate that 

r]t{x) = Hm-ix) >e} + Bn{x) l{rit-{x) = 9}. 

In the proof of both theorems, we first study the process {£,t) that includes mixing events but 
excludes local events whose dynamics is therefore described by the Markov generator 

hfio = E (5) 

Note that this process can be constructed graphically as previously by only using the Poisson 
processes attached to the edges of the graph. Note also that, since the state at each vertex is a 
convex combination of the states of the vertices at earlier times, 

^s(x) > 6 for all X G y implies that £,t[x) > 9 for all {x,t) £Vx (s,oo) 

and the analogous implication obtained by flipping the inequalities. Following the terminology 
introduced in 0], we call respectively upper/lower configurations the sets 

n+ = {^-.V ^ [0, 1] such that ^(x) > 6 for ah x e V} 
n_ = {^-.V [0, 1] such that ^(x) < 9 for all x e V} 

and observe that the previous implication means that, once the process ([5]) hits the set of upper 
configurations, it stays in this set forever. By definition of the Allee threshold, the same holds for 
the original process ([T]) from which we deduce the following lemma. 

Lemma 3 - For the process ^ on a finite graph, 

rjt = 1 for some t if and only if rjt € 1^+ for some t 
r]t = for some t if and only if r]t € r2_ for some t. 

Lemma [3] is one of the keys to proving both theorems. According to the lemma, it suffices to prove 
that, with the appropriate probability, the process on the complete graph hits a lower configuration 
whereas the process on the ring hits an upper configuration. 

4. The process on the complete graph 

This section is devoted to the proof of Theorem [TJ As previously mentioned, the first step is to 
study the process ([5]) that excludes local events but includes mixing events. It is obvious that, when 
starting with a single vertex in state 1, this process eventually hits a lower configuration provided 
the number of vertices is sufficiently large. The key to the proof is to show that the time to hit 
a lower configuration can be made arbitrarily small, which relies on large deviation estimates for 
several random variables that we now define. We call collision a mixing event that involve two 
occupied vertices, i.e., a collision occurs at time t whenever 



t = Tn{x,y) and min(^t_(x), / 0. 
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Then, we define the three random variables 

time to dispersion : td := inf {t : £ 
time to colhsion : tq first time a colhsion occurs 
number of occupied patches : := card {x G Z/A^Z : ^((x) 7^ 0}. 

We will prove that, with probability close to one when is large, 

TD < Tn and Tc > Tn and < 4^^^^ (6) 

for some T/v that tends to zero as — >■ 00. The probability of the events in ([6]) will be estimated 
backwards by conditioning, i.e., the probability involving the time to collision is obtained by con- 
ditioning on the number of occupied patches, while the probability involving the time to dispersion 
is obtained by conditioning on the time to collision. To complete the proof, we will return to the 
process ([1]) and use our estimates for the probability of the first and last events in ^ to prove 
that the probability that a local event occurs in any of the occupied patches before the time to 
dispersion tends to zero as the number of vertices — )• 00. 

Mapping to a dynamic graph - To estimate the probability of the events in (l6|), we first define 
a mapping to visualize the evolution of the process ([5]) through a dynamic graph, i.e., a continuous- 
time Markov chain whose state at time t is a random oriented graph 

Ht := {Vt,Et) where Vt C F x N. 

The dynamic graph is coupled with the process ([5]) and defined as follows. 

• The graph Hq has only one vertex, namely (xq ,0) G V x N where xq is the single vertex in 
state 1 initially, and no (oriented) edge. 

• We can ( leaf at time t whenever (x, i) G Vt and (x, i + 1) ^Vt- 

• Assume that (x,z) G Vt- is a leaf and that t = Tn{x,x') for some n > 1. Then, we define the 
new vertex set and the new edge set as 

Vt := Vt- U {(x,i + l),(x',i + l)} 

Et := Et- U {{x,i) ^ {x,i + l),{x,i) ^ {x',i + l)}. 

The left-hand side of Figure [2] gives a schematic picture of a realization of the process with only 
mixing events where Xj denotes the jth vertex that becomes occupied, while the right-hand side of 
the figure gives a picture of the corresponding graph. Note that a vertex has a positive density at 
time t if and only if it is the first coordinate of a leaf of Ht- It follows that 

l^tl = number of leaves in Ht (8) 

and, using in addition the evolution rule ([7]), that a collision event results in two pairs of oriented 
edges each pointing to the same leaf. In particular, 

Ht is an oriented binary tree with root (xq, 0) if and only if t < rc. (9) 

From properties ([7]) and ([9]), it also follows that 



{x,i) is a leaf at time t < tc implies that (,t{x) < (1 — ^)*. 



(10) 




Figure 2. Pictures related to the proof of Theorem [T] 



In fact, if (x,i) is a leaf at time t < tq then there exists a unique oriented path from the root to 
this leaf and the density at vertex x can be computed explicitly looking at the number of vertical 
edges in this path, but this property is not needed in the proof of the theorem. 

Time to dispersion - The next step is to use ([8|)- (jl0p and the dynamic graph representation of 
the process to estimate the probability of the three events in Under the assumptions of the 
theorem, there exists n such that (1 — < 9. We then define 

Tn := n In {In N)/N and Kn := 4"''^(''^^) = 4^^^. 

The following three lemmas give estimates of the probability of the last event, the second event 
and the first event in ([6|), respectively, for the deterministic time T/v defined above. 

Lemma 4 - There is a > such that P {\Ctjs, 1 > I^n) < (In iV)^" for aU N large. 

Proof. The number of leaves is maximal when there is no collision, in which case the number of 
leaves jumps from i to i + 1 at rate iN. This together with ([5]) implies that 

E\Ctn\ = (number of leaves in i?Tjv) < 2^^^ = 
In particular, large deviation estimates for the Poisson distribution give 

P{\^tJ>Kn) = Pd^T^j >4^^^) < exp(-aln(lniV)) = (IniV)-'^ 
for a suitable constant a > and all sufficiently large. □ 
Lemma 5 - For aU N large, P {tq <Tn) <2 (IniV)"'*. 

Proof. Given that the graph Ht has i leaves, the probability of a collision at the next update of 
the system is equal to i/N . In particular, the conditional probability of a collision before T/v given 
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that the number of leaves at that time is smaller than Kn is 



N - 2N 

i=l 



< 



exp(4nln(lniV)) (IniV) 



4n 



2N 2N 



This, together with Lemma [H implies that 

P{tc<Tn) < P{rc<TN\\CTj<KN) + P{\^tJ> Kn) 
< (ln7V)4"/2iV + (IniV)-'* < 2(lniV)-'^ 

for all A'" sufficiently large. □ 

Lemma 6 - Let o > as in Lemmas [H and [5l Then, for all N large, 

P{td>Tn) < 2"+i(lniV)-^ + 2(lniV)^''. 

Proof. Motivated by pil|) . we introduce the stopping times 

(Tj := inf {t : for each leaf (x, i) € Vt we have i > j}. 

The first step is to prove that (T„ < T/v with probability arbitrarily close to one when the number 
of vertices is large. Note that, according to the evolution rules ([7]), we have 

card {xeV : (x, i) G Vt} < 2' for all t > 0. (11) 

Moreover, since each vertex has degree — 1 and is therefore involved in a mixing event at the 
arrival times of a Poisson process with intensity — 1, we have 

P{mf{t: {x,i + 1) £Vt}- inf {t:{x,i) eVt}>T) = exp(-(iV - 1) T) (12) 

for all X € y such that (x,i) € Vt at some time t, i.e., the amount of time a vertex in the 
dynamic graph is a leaf is exponential with parameter — 1. From (jllh and ()12p . we deduce that 
the temporal increment — aj required to grow one more generation in the dynamic graph is 
stochastically smaller than the maximum of 2^ independent exponential random variables with the 
same parameter A^ — 1. In particular, having a collection ei,e2,e3, ... of independent exponential 
random variables with parameter A^ — 1, we deduce that 

n-l 

P {an >Tn) < P {aj+i - a J > Tn/tl) 

n-l n-l 

< ^ P (max {ei:i< 2^} >TN/n) < P (e* > T^/n for some i < 2^) (13) 
i=o j=o 

< ^ 2^' exp(-(A^- 1) ln(lnAf)/Af) < 2"+^(lnA^)-^ 
i=o 

for all A^ large. In other respects, using (llOh and recalling the definition of n, 

Cn < < Tc implies that Ctn{^) < (1 — Z^)" < ^ for each leaf {x,i) € Vtjv 

(14) 

implies that to < Lat. 
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Combining (I13p and (I14p . and using Lemma [SI we conclude that 

P{td>Tn) < P{aN>TN) + P{rc<TN) < 2"+i(lniV)-i + 2(lniV)-'^ 

which completes the proof. □ 

To conclude the proof of Theorem [H we return to the process ([1]). The proof is based on the 
simple observation that processes ([1]) and ([5]) are equal as long as no local event occurs in any of 
the vertices not in state 0. The bound on the number of occupied patches and the bound on the 
time to dispersion given respectively by Lemmas U] and [6] show that the time to dispersion for the 
process that includes local events is small as well with probability close to one, so the result follows 
from Lemma [3j This argument is made rigorous in the following lemma. 

Lemma 7 ~ For all N sufficiently large, we have 

pU0,Ij) < 2"+i(lniV)-i + SilnN)-" + n (lniV)2"+ViV. (15) 

Proof. The two processes ([TJ and ([5]) being constructed from the same graphical representation 
are equal as long as no local event occurs in any of the vertices not in state 0. Since local events 
occur at each vertex at rate one, having a collection ("1, C2, Cs, • • • of independent exponential random 
variables with parameter one, we deduce that 

P {r]t ^ Ct for some t <Tn\ \^Tn I < Kn) 

< P{mm{Q ■.i<KN}< Tn) = 1 -P(min{Ci : i < Kn} > Tn) (16) 

< 1 - exp (-i^TvTAr) < 1 -exp (-n (In 7V)2" ln(ln iV)/iV) < n (In iV)2"+i/iV 

for all N large. From Lemmas H] and [6l and ()16p . we obtain 

P ir]T^ in^) < P ^n^) + P {7]t ^ it for some t < Tn) 

< P{td> Tn) + P{r]t^ 6 for some t < Tjv [ [Cr^vl < Kn) + P {\iTj > Kn) 
< 2"+i(lniV)-i + 3(lniV)-'^ + n (In iV)2"+i/iV 

for all large. Since according to Lemma [3] we have 

p+{e, fi) = P {r]t G VL+ for some t) = P {rjt ^ ft^ for all t) < P {r]Tj^ ^ n^) 

the proof is complete. □ 

The theorem directly follows from Lemma [7] by observing that the right-hand side of (jlSp tends to 
zero as the number of vertices goes to infinity. 

5. The process on the ring 

This section is devoted to the proof of Theorem [2l To understand the process on the ring, the 
first step is to study its counterpart on the infinite one-dimensional lattice using a so-called block 
construction. The idea is to couple a certain collection of good events related to the infinite system 
with the set of open sites of a one dependent oriented site percolation process on 



7i := {{x, n) € Z X : x + n is even}. 
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For a precise definition and a review of oriented site percolation in two dimensions, we refer the 
reader to Durrett [Ij. This couphng together with results from [2j implies that, starting with a single 
occupied patch, there exists with positive probability a linearly expanding region that contains a 
positive density of patches above the Allee threshold. The second key step is to prove that in fact 
all the patches in this space-time region are above the Allee threshold, from which it follows that, 
with the same positive probability, the process on the finite ring starting with a single occupied 
patch hits an upper configuration before it hits a lower configuration. To prove linear expansion in 
space of the set of patches that exceed the Allee threshold, we observe that, under the assumptions 
of the theorem, there exists a constant a fixed from now on such that 

a < (1-//)^^ and a fi"^ {I - fif^ > 6 where T := 95. 

To define our collection of good events, we also introduce random variables that keep track of the 
number of mixing events and local events in certain space-time regions of the graphical represen- 
tation. More precisely, we introduce the number of mixing events 





:= card{n 


Tn{3,j + l) 


G(0,r)} fori = -1,-2 




:= card {n 


Tn{j-lJ) 


G(0,T)} forj = +l,+2 




:= card {n 


TnU-lJ) 


G {T,2T)} 



+ card {n : T„(i, j + 1) G (T, 2T)} for j = -1, +1 

as well as the number of local events 

Zj := card {n : G (T, 2T)} for j = -1, +1. 

From these random variables, we define the good event 

Q := {min(A:_i,A:i,Z_i,Zi) / 0} n {max(X_2,^2) < 2T} n {max(y_i,yi) < AT} 

as depicted in Figure [31 For every (x,n) G 7i, we define the good event Q{x,n) similarly but from 
the graphical representation of the process in the space-time region 

R{x,n) := {x -2,x + 2) X {2nT,2nT + 2T). 

The motivation for introducing these events is that, conditioned on Q{x, n), if the population density 
at patch X at time 2nT exceeds a then the same holds for the two adjacent patches 2T units of time 
later. By translation invariance of the evolution rules of the process in space and time, it suffices 
to prove the result for x = n = 0, which is done in the following lemma. 

Lemma 8 - We have {??o(0) > a} n 0(0,0) C {min (r/2T(-l), ?/2t(1)) > a}- 

Proof. The first step is to prove that the event on the left-hand side is included in the event that 
the population density at —1 and 1 exceed the Allee threshold at time 2T for the process ([5]) that 
includes mixing events but excludes local events. Assume that 

s_i := ri(-l,0) < ri(0,l) =: Si < T. 

Note that the last inequality si < T follows from Xi ^ 0. To study the process up until time si, 
we introduce the number of mixing events 



M := card{n : r„(-2,-l) G (s„i,si)}, 
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Figure 3. Picture of the good event fi = fi(0, 0). 



and in case M 7^ 0, the time of the first mixing event 

s := inf({T„(-2,-l) :n> 1} n 

Since the population density at patch —1 and at patch between time s_i and time s are convex 
combinations of their counterpart at time s_i, we have 

min(^t(-l),^t(0)) > min(^s_i(-l),^s_,(0)) > min (a/x, a(l - /x)) = a/x 
for all t G (s-i, s). In particular, 

min(e^(-l),^^(0)) > afi{l-fi). 
Since M < X_2 < 2T, we deduce from a simple induction that 

min(6(-l),6(0)) > a^^{l-^i)^ > a/x(l-/x)2^ > 9 
for all t G (s-i, si), and by definition of si, 

61 (1) = /x6i-(0) + (l-/x)6i-(l) > /^6i-(0) > a/x2(l-/x)^. 
In particular, using again our inductive reasoning and the fact that 

card{n:r„(-2,-l) G (si,T)} 

+ card{n : r„(l,2) G (si,T)} < X_2 + ^2 - M < 4T - M 

we deduce that 

min(et(-l),et(0),6(l)) > (l-/x)4^-^ min (61 (-1), 61(0), 61(1)) 



(17) 



(18) 



(19) 
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for all te{si,T). Finally, since Y-i < 4T and Yi < AT, 

min(6(-l),6(l)) > (l-Z^r min(^T(-l),eT(l)) 

> (1-^)4^ a/i2(l-^)4T = a/i2(l-^)8r > 9 

for all t G (T, 2T). Returning to the system with local events, since 

Ct(0) > e for ah t e (0, T) and ^f(±l) > 9 for ah t G (s±i, 2r) 

according to pT|) ~ p^ . these inequalities remain true for the original process ([T]). In particular, the 
populations at patches —1 and 1 exceed 9 between times T and 2T therefore 

r]t{—l) = 1 for some t E {T,2T) and = 1 for some t € {T,2T) 

since Z_i / and Zi / 0. Using again that Y-i,Yi < 4T, we conclude that 

mri-l) > > a and 7?2t(1) > (1 - m)^^ > a 

which completes the proof of the lemma. □ 

To deduce from Lemma [8] the existence of a linearly expanding region with a positive density 
of patches above the Allee threshold, we now prove that the common probability of all our good 
events exceeds the critical value pc of one dependent oriented site percolation. 

Lemma 9 - For T = 95, we have P n)) > 1 - 3^36 > p^. 
Proof. The key is simply to observe that 

• X^2^X^i,Xi and X2 are Poisson random variables with parameter T, 

• y_i and Yi are Poisson random variables with parameter 2T and 

• Z_i and Zi are Poisson random variables with parameter T. 

Using in addition that these random variables are independent, we deduce that 

P{9.{x,n)) = P{n) > l-AP{Xi = 0)-2P{X2>2T)-2P{Yi> AT) 

= 1 _ 4e-T _ 2 y e-^ _2 y l^^L g-^^ > 1 - S^^e 

n! n! 

n>2T n>AT 

when T = 95. The second inequality pc < 1 — 3^^^ in the statement follows from the contour 
argument described in, e.g., [1], Section 10. This completes the proof of the lemma. □ 

Lemmas [8] and [9] and the fact that 

R{x,n) n R{x',n') = whenever \x — x'\ > 1 or n ^ n' 

are the assumptions of Theorem 4.3 in Durrett pj with M = 1 and 7 = 3~^^ from which it follows 
that, for the infinite system starting with a single occupied patch at the origin, 

Wn '■= {x £ Z : (x, n) gV. and ri2nT{x) > a} 
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dominates stochastically the set of wet sites Wn at level n of a one dependent oriented site percola- 
tion process with parameter 1 — 7 and initial condition Wq = Wq. Since 1 — 7 > Pc we deduce that, 
with positive probability at least equal to the percolation probability, the set of patches that exceed 
the constant a expands linearly. This only proves persistence of the metapopulation of the infinite 
lattice, which is not sufficient to deduce global expansion of the system on the ring. The last step 
is to show that all patches in the expanding region exceed the Allee threshold. More precisely, on 
the event that percolation occurs, i.e., Wn / for all n, we have 

lim In '■= lim min Wn = —00 and lim rn '■= lim max Wn = +00 (20) 

n— >-oo n— ^-oo n— >oo n^oo 

and thinking of the infinite system as being coupled with one dependent oriented site percolation 
in such a way that Wn C Wn for all n, we have the following lemma. 

Lemma 10 - Assume that Wn ^ for all n. Then, for all 2; S Z, 

r]t{x) > 9 for all times t sufEciently large. 

Proof. Recall from the proof of Theorem 4.3 in Durrett |2] that the processes are coupled in such 
a way that the set of open sites for the percolation process is included in the set of good sites, 
where site (x,n) is said to be good whenever the good event Cl{x,n) occurs. This, together with 
the definition of the right edge, implies that there is a good path from site (0, 0) to site (r^, n), i.e., 
a sequence of integers xo = 0, xi, . . . , Xn = Tn such that 

(xm, "i) is good for m = 0, 1, . . . , n and \xm — Xm-i\ = 1 ior m = 1,2, . . . ,n. 

Since 7/o(0) > a, it follows from (|17p -(|19 p in the proof of Lemma[8]that, for m = 0,1, . . . ,n. 



rit{xm) > 9 for ah t e [2m,T,2mT + T) 
r]t{xni + l) > for ah t e [2mT + T,2mT + 2T). 

Similarly, there is a path (xq, 0) — ?> (x_i, 1) ^ • • • — > n) = n) such that 

r]t{x-rn) > for ah t e [2mT,2mT + T) 
r]t{x^„, -1) > e for ah t G [2m,T + T, 2mT + 2T) 



(21) 



(22) 



for m = 0,1, . . . ,n, and we may assume that x-m < Xm for all m. We claim that all patches in the 
space-time region delimited by (f2T]) and ([22]) are above the Allee threshold 9, i.e., 

r]t{x) > for all {x,t) G [x^m,Xm] x [2mT,2mT + T) 

r]t{x) > e for aU (x, t) G [x_m - 1, x^ + 1] x [2mT + T, 2mT + 2T), ^'^'^^ 



which we prove by induction. Assume that (I23p holds for some m < n. The fact that this again 
holds at time t = 2mT + 2T simply follows from the fact that 

[x_m_i,Xm+i] C [x_m-l,Xm + l] since X_m-1 > - 1 and Xm+l<Xni + l- 



To prove that this holds at later times, we distinguish three types of events: 
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• Local events cannot violate the first line of ()23p since patches above the Allee threshold can 
only experience local expansions to their carrying capacity. 

• Mixing events in [x-m-i,Xm+i] cannot violate the first line of (j23p since the new states of 
interacting patches above the Allee threshold are again above the Allee threshold. 

• Mixing events along i 1)2^— m—i 

) or {xm+i,Xm+i + 1) can violate the first line of ([23]) 
but this would contradict either ([2T]) or ([22]) . 

This proves the first line of (|23p at step m + 1 while the second line follows from the exact same 
reasoning. From ([23]) . we deduce that 

r]t{x) > e for all (x, t) G x [2mT, 2mT + 2r). 

In particular, the lemma follows from ([20p . □ 

Returning to the process on the ring, we deduce from Lemma [TU] that, with positive probabil- 
ity at least equal to the percolation probability and starting with a single patch in state 1, the 
system reaches an upper configuration before it reaches a lower configuration, an event that leads 
to global expansion according to Lemma [3l This completes the proof of Theorem [2j 
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